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Abstract 

A search for pair-produced doubly charged Higgs bosons has been performed 
using the data collected by the DELPHI detector at LEP at centre-of-mass 
energies between 189 and 209 GeV. No excess is observed in the data with 
respect to the Standard Model background. A lower limit for the mass of 97.3 
GeV/c 2 at the 95% confidence level has been set for doubly charged Higgs 
bosons in left-right symmetric models for any value of the Yukawa coupling 
between the Higgs bosons and the r leptons. 
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1 Introduction 

Doubly charged Higgs bosons (i? ±± ) appear in several extensions to the Standard 
Model [1], such as left-right symmetric models, and can be relatively light. In Super- 
symmetric left-right models usually the SU(2)r gauge symmetry is broken by two triplet 
Higgs fields, so-called left and right handed. Pair-production of doubly charged Higgs 
bosons is expected to occur mainly via s-channel exchange of a photon or a Z boson. In 
left-right symmetric models the cross-section of e + e~ — > H£ + H£~ is different from that 
for e + e~ — > H^ + H R ^, where Hf^ and are the left-handed and right-handed Higgs 
bosons. The formulae for the decays and the production of these particles can be found 
in [2]. 

In these models the doubly charged Higgs boson couples only to charged lepton pairs, 
other Higgs bosons, and gauge bosons, at the tree level. The current limit and the mass 
range of this analysis is restricted to the interval between 45 GeV/c 2 , the LEP1 limit 
set by OPAL [3], and the kinematic limit at LEP2, that is around 104 GeV/c 2 . The 
dominant decay mode of the doubly charged Higgs boson is expected to be a same sign 
charged lepton pair, the decay proceeding via a lepton number violating coupling. As 
discussed in [2], due to limits that exist for the couplings of i/ ±=t — > e ± e ± from high 
energy Bhabha scattering, i7 ±± — > /x ± jU ± from the absence of muonium to anti-muonium 
transitions and H ±i: — * / u ± e ± from limits on the flavour changing decay /i ± — > e T e ± e ± , 
electron and muon decays are not likely. In addition, most of the models expect that the 
coupling to rr will be much larger than any of the others. Therefore, only the doubly 
charged Higgs boson decay if ±± — > r ± r ± is considered here. 

The partial width for the H ±i: decay into two r leptons is, at the tree level [2]: 

r ( H ±±^+±\ Kt^ 1 1 2m ?Wi 4m 2 . l /2 

where m T is the mass of the r lepton and h TT is the unknown Yukawa coupling constant. 
Depending on the h TT coupling and the Higgs boson mass the experimental signature is 
different. If h TT is sufficiently large, h TT > 1CT 7 , the Higgs boson decays very close to the 
interaction point. We describe here an analysis to search for such events. If h TT is smaller 
the decay occurs inside the tracking detectors or even beyond them, making this analysis 
inefficient. In this case pre-existing analyses were applied which are further discussed 
below. 



2 Data sample and event generators 

The data collected by DELPHI during the LEP runs at centre-of-mass energies from 
189 GeV to 209 GeV were used. The total integrated luminosity of these data samples 
is ~ 570 pb" 1 . The DELPHI detector and its performance have already been described 
in detail elsewhere [4,5]. 

Signal samples were simulated using the PYTHIA generator [6]. In this analysis 
samples with doubly charged Higgs boson with masses between 50 and 100 GeV/c 2 , 
in 10 GeV/c 2 steps, were used at different centre-of-mass energies, both for left-handed 
and right-handed bosons, and different Yukawa coupling constants. 

The background estimates from the different Standard Model processes were based on 
the following event generators, interfaced with the full DELPHI simulation program [5]. 
The WPHACT [7] generator was used to produce four fermion Monte Carlo simulation 
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events. The four fermion samples were complemented with dedicated two photon colli- 
sion samples generated with BDK, BDKRC [8] and PYTHIA [6]. Samples of qq("j) and 
// + //~(7) events were simulated with the KK2f generator [9]. Finally, KORALZ [10] was 
used to simulate t + t~(^) events and the generator BHWIDE [11] was used for e + e~{^) 
events. 

3 Data selection 

The search for pair-produced doubly charged Higgs bosons makes use of three different 
analyses depending on the h TT coupling or, equivalently, on the mean decay length of the 
Higgs bosons. When the mean decay length of the Higgs boson is very small, the resulting 
final state consists of four narrow and low multiplicity jets coming from the interaction 
point. This analysis is explained in detail in section 3.1. For intermediate mean decay 
lengths of the Higgs boson the topology consists of two tracks coming from the interaction 
point, and with either secondary vertices or kinked tracks. If the Higgs boson decays 
outside the tracking devices the signature corresponds to stable heavy massive particles. 
These two analyses were designed for the search for supersymmetric particles decaying 
to similar topologies. Details can be found in [12]. 

3.1 Small impact parameter search 

An initial set of cuts was applied to select events with four jets of low multiplicity. 
Only tracks with an impact parameter below 4 cm both in the plane transverse to the 
beam axis and in the direction along the beam axis were considered in the analysis. A 
charged particle multiplicity between 4 and 8 was required. Events were clustered into 
jets using the LUCLUS algorithm [6], requiring each jet to be separated from the others 
by at least 15 degrees, and only events with four reconstructed jets were accepted. To 
improve the reconstruction of the r energy, the r momenta were rescaled, imposing energy 
and momentum conservation and keeping the r directions at their measured values. If 
the rescaled momentum of any jet was negative, the event was rejected, as such events 
are commonly not genuine four jet events. 

The two photon background was reduced by the following energy and momentum 
requirements: the energy of observed particles produced at a half opening angle to the 
beam axis exceeding 25° had to be greater than 0.15\/i, the momenta of the jets were 
required to be larger than O.Oly^ and the total neutral energy had to be less than 0.35a/s. 

The four lepton background was rejected by requiring that the momentum of the most 
energetic lepton identified (electron or muon) was less than 0.25^ and the momentum 
of the second most energetic lepton identified was less than 0.15\/s. The algorithms used 
in the lepton identification were the same as those used in the selection of fully-leptonic 
W pairs [13]. 

The calculated r momenta, defined above, were used to reconstruct the Higgs boson 
mass. The charge of the r jet was calculated as the sum of the charges of its constituent 
particles. If this value was not ±1, then the charge of the most energetic charged particle 
was assumed to be the charge of the r. For events with two positive r lepton candidates 
and two negative r lepton candidates the charge was used to assign the pairing of both 
doubly charged Higgs bosons. If the total charge was not equal to 0, the pairing was chosen 
to minimise the difference between the two reconstructed masses of the Higgs bosons. The 
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ratio , ' , Rec++ Rec — r-L W as required to be less than 0.7. Finally the reconstructed event 

( M Rec + ++ M Rec—)/ 2 ^ J 

mass, defined as the average of the two masses, had to be greater than 40 GeV/c 2 . 

The effects of the selection cuts are shown in Table 1 for the combined 189-209 GeV 
sample. After all cuts were applied only one event was observed in the data with a mass of 
69±3 GeV/c 2 , while 0.9 events were expected from background processes. The candidate 
was collected at v /i=206. 7 GeV and is compatible with the assignment ZZ — > t + t~t + t~ . 
The most probable reconstructed masses with different sign leptons are indeed compatible 
with a Mz-Mz mass hypothesis at the one sigma level. The signal efficiency was around 
40% for a wide range of masses between 70 and 100 GeV/c 2 for both left-handed and 
right-handed doubly charged Higgs bosons, as shown in Table 2. Table 3 shows the 
selection efficiencies for left-handed doubly charged Higgs bosons for several if ±± masses 
and several h TT couplings at ^=206.7 GeV. The final reconstructed mass spectrum and 
the expected mass distribution in simulated signal events are shown in Figure 1. The 
good level of agreement between data and simulation observed at different stages of the 
analysis is demonstrated in Figure 2. 



cut 


data 


total bkg. 


an 


other 


£ H++H-- 


Four jets preselection 


59 


67.41±0.95 


44.01±0.31 


23.40±0.90 


59.2% 


anti 77 cuts 


26 


31.03±0.48 


28.90±0.25 


2.13±0.41 


52.3% 


anti 4 lepton cuts 


1 


1.87±0.07 


1.69±0.06 


0.18±0.03 


48.7% 


Mass requirements 


1 


0.91±0.04 


0.85±0.04 


0.06±0.01 


44.2% 



Table 1: The total number of events observed and the expected background after the 
different cuts used in the analysis for the small impact parameter search for the com- 
bined 189-209 GeV sample. The errors are only statistical. The last column shows the 
efficiency for a left-handed doubly charged Higgs boson signal with m H ±± = 100 GeV/c 2 

at ^=206.7 GeV. The statistical error in the signal efficiency is about 1.5% in all cases. 



channel 


M H ±± (GeV/c 2 ) 




50 


60 


70 


80 


90 


100 


left-handed 


32.7 


36.6 


40.5 


44.8 


43.4 


44.2 


right-handed 


31.8 


37.0 


40.0 


44.0 


44.8 


45.2 



Table 2: Selection efficiencies (in %) for left-handed and right-handed H ++ H — > 
t + t + t~t~ for several H ±i: masses and h TT > 10~ 7 at ^=206.7 GeV, for the small 
impact parameter search. The statistical error is around 1.5% in all cases. 



3.1.1 Systematic uncertainties 

Several sources of systematic uncertainties on the signal efficiency and the background 
level were investigated. The particle identification was checked on di-lepton samples both 
at the Z peak and at high energy. The discrepancy in the efficiencies between the data 
and the simulation was found to be lower than 2% in all cases. The track selection and the 
track reconstruction efficiency was also studied with these samples. These effects were 
studied by the comparison between data and simulation for tracks at the boundaries 
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h TT 


M H ±± (GeV/c 2 ) 


50 


70 


90 


100 


4- 10- 8 


0.2/38.1/13.1 


1.6/43.0/1.4 


6.0/23.9/0.0 


20.5/5.3/0.0 


io-« 


0.0/6.4/68.4 


0.0/16.0/57.2 


0.0/30.5/22.7 


0.0/36.3/7.3 


< 10- 9 


0.0/0.0/77.6 


0.0/0.0/77.6 


0.0/0.0/41.3 


0.0/0.0/41.6 



Table 3: Selection efficiencies (in %) for left-handed doubly charged Higgs bosons for 
several if ±± masses and several h TT couplings at ^=206.7 GeV, for the three analyses 
performed (small impact parameter search, search for secondary vertices or kinks and 
search for stable massive particles, respectively). The statistical error is around 1.5% in 
all cases. 

of subdetector acceptances, where systematic effects are expected to be larger. The 
systematic error of these effects was about 1.5%. 

The errors on the background and signal rates from the modelling of the detector 
response were a few percent. Different variables at preselection level have been studied, 
with good agreement between data and simulation observed. The distributions in relevant 
variables before the anti 77 cuts and the anti four lepton cuts are shown in Figure 2. The 
masses reconstructed from both same sign and different sign lepton pairs, before the anti 
four lepton cuts were applied, are shown in Figure 3. For the opposite sign lepton pairs 
only the mass of the combination closest to the Z mass has been given and the Z peak is 
visible. 

The total systematic error on the background was about 13%, with a dominant con- 
tribution of about 12% due to the limited simulation statistics available. The total 
systematic error on the efficiency was about 5%. 

3.2 Search for secondary vertices or kinks 

When the lifetime is such that the particle decays inside the tracking detector, the 
previous analysis is inefficient, because impact parameter cuts are applied to reject the 
background coming from secondary interactions. We have applied here the analysis de- 
scribed in [12], that performs a special track reconstruction for this particular topology, 
looking for decay vertices far from the interaction point. 

After all cuts five events were selected in the data, while 2.9 events were expected 
from the background. The signal efficiency was about 40%, if the mean decay length 
was about 50 cm with a smooth fall for both lower and higher mean decay lengths. The 
selection efficiencies for several if ±± masses and several h TT couplings at ^=206.7 GeV 
are shown in Table 3. 

3.3 Search for stable massive particles 

If the lifetime is even larger, the H ±:k crosses the tracking devices without decaying. 
The analysis described in [12] to search for stable heavy particles is applied here. It is 
based on the measurement of anomalous ionisation loss measured in the Time Projection 
Chamber and of the absence of Cherenkov light detected in the Ring Imaging Cherenkov 
Detector. 
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One event was selected in the data, in agreement with the expected background of 1.9 
events. For stable particle masses in the range of 50-80 GeV/c 2 the efficiency was ~ 75%, 
decreasing to ~40% for masses near the kinematic limit (Table 3). 

4 Determination of the mass limit 

No evidence for H ++ H production was observed. A modified frequentist likelihood 
ratio method [14] has been used to compute the cross-section and mass limits. The re- 
constructed event mass was used as a discriminant variable in the computation of the 
confidence levels in the small impact parameter analysis, while for the others only the 
number of events were used. The systematic errors were taken into account in the com- 
putation. All centre-of-mass energies and the three analyses were treated as independent 
experiments. For intermediate mean decay lengths of the Higgs bosons in many cases 
two analyses have significant efficiency. However the overlap of the samples selected by 
the analyses, both for the signal and for the background, was negligible. 

A very similar behaviour, both in terms of efficiency and of mass distributions, was 
observed for the left-handed and the right-handed doubly charged Higgs bosons. Hence, 
the average of both contributions were used to calculate the confidence levels. The 
expected left-handed and right-handed cross-sections were calculated using the PYTHIA 
generator [6]. 

Previous searches for H ±zL pair production have already excluded M#±± < 45.6 
GeV/c 2 [3]. Therefore, this search was limited to masses greater than this value. The 
limits at 95% confidence level for different values of h TT are shown in table 4. Figure 4 
shows the 95% confidence level upper limits on the cross-section at y/s = 206.7 GeV for 
the production of H ++ H — > t + t + t~t~ for these values of h TT . The comparison of 
these limits with the expected cross-section for left-handed and right-handed 
pair production yields 95% confidence level lower limits on the mass of the if^ ± and 
bosons of 98.1 and 97.3 GeV/c 2 , respectively for any value of the h TT coupling. 

This search slightly improves previous searches for h TT > 10~ 7 [15], and in addition is 
extended to the whole range of the h TT coupling. 



h 


Left-handed 


Right-handed 




Observed 


Expected 


Observed 


Expected 


> io- y 


99.6 


99.6 


99.1 


99.1 


4 • 10" 8 


98.1 


98.4 


97.3 


97.6 


10"« 


99.0 


99.4 


98.4 


98.9 


< io- y 


99.6 


99.6 


99.3 


99.3 



Table 4: Median expected and observed if ±± mass limits at 95% C.L. in GeV/c 2 for 
different values of the h TT coupling. 



5 Conclusion 

A search for pair-produced doubly charged Higgs bosons decaying into r leptons was 
performed using the data collected by DELPHI at LEP at centre-of-mass energies from 
189 GeV to 209 GeV in R-parity conserving supersymmetric left-right symmetric models. 
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Three different analyses were applied to cover the whole range of the h TT coupling: decays 
very close to the interaction point, inside the tracking detectors or beyond them. No 
significant excess was observed and a lower limit on the doubly charged Higgs boson mass 
of 97.3 GeV/c 2 has been set at 95% confidence level for any value of the h TT coupling. 
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Figure 1: The reconstructed mass distribution after all cuts for the small impact param- 
eter search. The hatched histogram corresponds to the expected background and the 
dot with the error bar shows the one remaining candidate event. The dotted line corre- 
sponds to simulated events with m u ±± = 70 GeV/c 2 and the dashed line corresponds to 

simulated events with m H ±± = 100 GeV/c 2 . 
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Figure 2: Event selection variable distributions at different stages of the analysis for 
the small impact parameter search. The top plots show the minimum momentum of 
the jets and the visible energy outside 25° around the beam axis scaled by ^/s after the 
four jet preselection cuts. The bottom plots show the momentum of the most energetic 
identified lepton and the momentum of the second most energetic identified lepton scaled 
by v^i after the anti 77 cuts. The solid lines show the expected background, the dots 
the observed data and the dashed lines correspond to m H ±± = 100 GeV/c 2 . The signal 
is multiplied by a factor 35 in the top plots and by a factor 4 in the bottom plots. 
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Figure 3: Reconstructed mass distributions for the small impact parameter search. The 
masses are shown for the same sign lepton pairs (top) and the opposite sign lepton 
combination closest to the Z mass (bottom). These distributions are shown before the 
anti four lepton cuts. The solid lines show the expected background, the dots the observed 
data and the dashed lines correspond to simulated events with m H ±± = 100 GeV/c 2 . 
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Figure 4: Upper limits, at 95% confidence level, on the production cross-section for a pair 
of doubly charged Higgs bosons as a function of the doubly charged Higgs boson mass 
at y / i=206. 7 GeV, assuming 100% branching ratio for the decay of if ±± into r ± r ± for 
different values of the h TT coupling. The dashed grey curve shows the expected upper 
limit with one and two standard deviation bands and the solid grey curve is the observed 
upper limit of the cross-section (the grey curves are those inside the bands). The dashed 
black and solid black curves show the expected production cross-section of and 
pairs in left-right symmetric models. 



